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Abstract 

This paper presents a mathematical model of a large class of machines which experience internal friction 
losses and maintain an essentially linear relation between input and output forces. The model is described in 
a general setting and a universal theorem is derived quantifying net work done and energy lost over a given 
process. Examples are given to illustrate the application of the model to predicting the performance of 
reciprocating heat engines equipped with such mechanisms. 

© 2003 Elsevier Ltd. All rights reserved. 


1. Introduction 

Force-linear mechanisms are machines in which the relation between input and output (or dual 
input) force or torque is as linear as it can be. Many real machines are exactly force-linear, ranging 
from simple harmonic speed reducers [1] to the crank-connecting rod mechanisms described be¬ 
low. Other machines can be sufficiently well approximated as force-linear, or as serial connections 
of such. This paper introduces a general mathematical model for force-linear machines and de¬ 
rives a formula for net work output in terms of a given input process. The formula classifies and 
quantifies how, when, and where these machines suffer frictional losses. 

The kind of insight this yields can be useful in designing and applying machines for specific 
tasks. To illustrate, the paper applies the force-linear results to reciprocating heat engines. It is 
shown first how the model can be used to evaluate the match between a mechanism and a 
thermodynamic cycle. A further application of the theory is the determination of the response of a 
large class of engines to pressurization or supercharging. Conditions are derived under which 
engine performance is improved by supercharging. At this point a counter-intuitive result is re¬ 
vealed and a specific example is constructed in which engine performance is actually degraded by 
supercharging. 
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2. Machines 

Fig. 1 schematically represents a basic 2-actuator machine. The machine, whose internals are 
not subject to analysis here, interacts with external entities through its actuators. Each actuator is 
a device capable of either linear or angular displacement, but not both, and capable of accepting 
or applying force or torque respectively. Coordinate systems, represented by x and y, are assigned 
to the range of motion of each actuator. The kinematic connection between the actuators is 
mathematically expressed with reference to a single parameter s. This is a parameter in terms of 
which the position of all of the machine elements can be specified. For example, s may be taken to 
be the crank angle in an analysis of a 2-stroke IC engine, or taken to be camshaft angle in a 4- 
stroke. In this paper, x(s) will denote the piston position and y(s) output shaft position when 
applied to engines. The operation or motion of a machine is represented by the time variation of 
the state parameter. 

2.1. Constrained motion 

For the purpose of analyzing machine operation in this paper, its motion in time is taken as 
fixed. The most common assumption for example in the analysis of an engine would be that its 
shaft has a uniform rotational motion. This is satisfied close enough if the flywheel is relatively 
large or there are a number of cylinders working on the same crankshaft. Motions other than 
uniform rotary can be handled as well, but the motion of the device undergoing analysis, whatever 
it is, is a priori specified. Thus an underlying assumption in any instance of analysis in this paper is 
that the parameter 5 changes in time in a prescribed way, that is, that a functional relation 5 = a{t) 
is given. 

The advantage of this constrained motion analysis is that it reduces the rest of the machine 
description to a relation solely between the actuator forces and the state parameter. With the 
motion in time of the machine state parameter specified, the motion of every part of the machine 
in time is set by its kinematic geometry. Taking the masses of the internal machine members into 
account, the force relation between / and g at a given state parameter value 5 results from a static 
type force analysis with joint friction and any appropriate inertia forces and torques included. 
Thus constrained motion analysis allows one to work with a force relation of the simple form 
F{s,f,g) = 0. 

It is important to note that if the operating regime s = a(t) changes, the force relation F(s,f, 
g) = 0 will change also. A given force relation remains in effect only as long as the motion giving 
rise to it remains the same. Note also that the machine state parameter can be taken to actually be 


y 



Fig. 1. Conceptual representation of a two-actuator machine. 
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clock time, s = t, but then one must be especially careful to remember that a fixed motion of the 
machine is assumed. A change of motion would entail a new parameterization of the machine. As 
long as this is kept in mind, it is actually simpler to use t as the state parameter [4], A machine with 
a prescribed motion will be denoted by M = [x(s),y(s),.F(,y,/,g)]. 

2.2. Sign conventions 

For uniformity, energy exchanges will be expressed relative to the mechanism. This means that 
although the orientation of the actuator scales can be arbitrarily chosen, once chosen, the positive 
external force or torque direction is taken to he in the direction of increasing actuator position co¬ 
ordinate. For simplicity it will be henceforth assumed in this paper that the motion of the 
mechanism is such that the state parameter 5 increases in time. There is no loss of applicability 
here, since a simple reparameterization can always be made to meet this requirement. 

Having o'(t) = s > 0 allows a straightforward identification of the direction of work flow. If/ 
represents the force or torque applied to the x-actuator, when / -x'{s) > 0, the force / is doing 
work on the machine at state s. If / • x'(s) < 0, then the machine is doing work against the external 
resistance / at state .v. 

2.3. Friction 

In its most basic role, the mechanism of an engine acts as a work transducer. It only exchanges 
energy with the rest of the engine in the form of mechanical work transmitted through the motion 
of its actuators and through the dissipation of some of this work in the form of heat of friction lost 
to the surroundings. 

2.4. Energy storage 

A secondary function performed by the mechanism of an engine is that of an energy store. This 
is an incidental capacity, arising primarily by virtue of the fact that machines have moving parts 
with mass, and so are capable of storing and releasing kinetic energy, Mechanisms can also 
involve compression of springs or the like, and so can also store and release potential energy 
internally. However, in typical engines, the internal energy content of the mechanism section is 
small compared to the energy stored in the flywheel, and the energy fluctuation of the internal 
members of an engine mechanism can usually be safely ignored. This state of affairs will be 
modeled here by postulating that the mechanism proper maintains a constant internal mechanical 
energy, that is the total kinetic and potential energy within the mechanism is constant throughout 
its operation. 

2.5. Law of the machine 

The first law of thermodynamics applied to the operation of a machine as modeled above re¬ 
quires that the work done by the machine never exceeds that done to it. The instantaneous ex¬ 
pression of this is: 
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Fig. 2. The force relation of a machine at a state s for which x'(s) and /(s) are both positive. 

if F(s,f,g ) = Oand.v > 0then fx\s) + gy'(s) ^ 0 (1) 

Fig. 2 depicts this condition. A hypothetical force relation F(s,f, g) = 0 is shown in the f-g plane 
for a particular value of s. Condition (1) is satisfied in as much as the graph of F is in one of the 
closed half-planes determined by the line g = -fx' (s)/y' (s). The graph is in the upper half-plane if 
y(s) is positive and in the lower if negative. Recall that s > 0 is assumed here; inequality (1) has 
the opposite sense when s < 0. 

The slope m v = x'(s)/y'(s) is called the velocity ratio at state .v. In a reciprocating engine, piston 
velocity relative to ,v, namely x'{s), alternates in sign while the output shaft typically has unidi¬ 
rectional velocity which can be taken to be positive. For such a machine, the force relation will 
always he above the velocity ratio line, g = —m v f. 

2.6. Force processes 

Flywheel equipped engines invariably have mechanisms with a special non-jamming charac¬ 
teristic. Whatever force (within a reasonable range anyway) the piston may apply to the mech¬ 
anism in any state, there must, in principle at least, be a corresponding torque which the flywheel 
could apply to sustain the operation of the engine. This means that for every state s and force 
/, there is a force g such that F(s,f,g) = 0. Typically g is uniquely determined and thus in 
applications to engines and most machines one may regard g as functionally dependent upon .v 
and /. 


3. Force-linear machines 

A machine operating with a given motion will be called force-linear or a v-machine if for each 
state s, the graph of F(s,f,g) = 0 in the f-g plane is the union of two coterminal rays. Fig. 3 
illustrates the force relation of a r-machine at a state where x'{s) and y'(s) are positive. 
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Fig. 3. The force relation of a r-machine where x'(s), /(s) > 0. 

Because of the restriction in this paper to non-jamming mechanisms, there is always a left and a 
right-hand ray for r-machines. The slopes of these rays will be denoted by M and M + respec¬ 
tively. These slopes are functions of the state s, and are typically piecewise continuous in s. The 
vertex ( f 0 ,g 0 ) = if o{s) , g o{s)) is also a piecewise continuous function of s. Note that for force- 
linear machines, the first law (1) can be written as 

if - faff) + (g~ go)y'{s) ^ 0 (2) 

where s > 0 is being assumed throughout as stated. 

A compact expression for a the force relation of a non-jamming force-linear machine is the 
following 

g - go = M + {f-f a ) + - M_{f -f o y (3) 

The positive and negative part functions used here are as in standard mathematical practice [2] 
namely: 

z + = max{0,z} and z~ = max{0, — z} 

Note that both parts are non-negative by definition and that the following is an identity: 

z = z + — z~. 

A common force-linear mechanism is the crank-connecting rod system shown in Fig. 4. With 
pure Coulomb friction on the main shaft, the crankpin, and the piston slider, the relation between 



Fig. 4. The crank-connecting rod mechanism. 
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Fig. 5. The force relation of a crank-connecting rod mechanism shown for every 30° of crank rotation. The horizontal 
axis represents the net force acting on the piston and the vertical axis represents the corresponding torque on the main 
shaft. Velocity ratio is represented by a dotted line. 

the force / acting on the piston and the torque g applied to the crankshaft is force-linear. Fig. 5 
shows the force relation of this mechanism for every 30° increment of crank rotation; large 
friction coefficients ji and large journal radii r were chosen in this example to more emphatically 
show the nature of the resulting force relation. The values of the parameters used in this example 



Fig. 6. Graphs of M + and M_ for the crank-connecting rod mechanism with friction coefficients set at 0.15. The central 
curve is the velocity ratio. With no friction, M + and M_ would coincide with the velocity ratio m v . 
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are: friction coefficients n s = i l P = R c = 10; connecting rod to crank throw ratio L/R = 4; and 
journal radius to crank throw ratios r c /R = 0.32 and r p /R = 0.28. 

For smaller friction coefficients, the graphs are similar but the /'-arms are closer to the velocity 
ratio line. If constant preloading is added to the piston or the main crankshaft, the relation re¬ 
mains force-linear but the vertex is moved away from the velocity ratio line as generically rep¬ 
resented in Fig. 3. Fig. 6 shows plots of the slopes M + and vs. crankangle s for the same 
mechanism but with the friction coefficients reduced to 0.15. 

In the following sections, a formula is derived for the shaft work output of a force-linear 
machine when subjected to a given force process. The form of the terms in the formula make it 
simple to identify and study the nature and magnitude of the friction losses incurred. 


4. Kinematic power ratio 

It is conceptually helpful to identify the ratio of net forces on the actuators normalized by the 
velocity ratio. This will be referred to as the kinematic power ratio. For force-linear machines there 
are two such ratios, denoted by k + and k , which correspond to the two rays of the /'-relation. The 
kinematic power ratios for a force-linear machine are defined by 

f -M a y/x' if x' y o , 

a \0 if x' = 0 1 j 

where a e {+,—}■ The kinematic power ratio is the ratio of the instantaneous power being done 
by the v-actuator to the power being delivered to the x-actuator. Note the minus sign in the 
definition. For applications to engines for example, if the x-actuator is the piston and the y- 
actuator is the output shaft, then the kinematic power ratios as defined by (4) represent the power 
out of the shaft divided by the power delivered to the piston. 

The ratios k + and k_ are functions of the state parameter s and can vary over the entire range 
of real numbers. Fig. 7 shows plots of k + and k_ for the crank-connecting rod mechanism of 
Fig. 6. Fig. 8 shows k + and k for the same mechanism with all the friction coefficients increased to 
1 . 0 . 



Fig. 7. Graphs of k + and for the crank-connecting rod mechanism of Fig. 6 with friction coefficients of 0.15. 
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Fig. 8. Graphs of k + and k _ for the crank-connecting rod mechanism as in Fig. 7 with the friction coefficients increased 
to 1.0. 

5. The work formula for force-linear machines 

Practical machines are evaluated by comparing the total energy exchanges between the two 
actuators of the machine performing a given task. Work is put into one end of the machine and 
the interest is in how much one gets out the other side. This is obviously the case with devices such 
as jacks, hoists, and automotive transmissions where effort is put into one side and beneficial effect 
appears out of the other. But it is also relevant when work alternately goes both out of and into 
each side as is the case with most reciprocating heat engines. Work is put into the piston side of 
the engine mechanism by the thermodynamic working substance, but also can be taken out at 
different times during the cycle to perform compression processes for example. The same is true on 
the flywheel side of the engine mechanism. Work is alternately put into and taken out of the 
flywheel as the engine operates. The net work out of the flywheel side over a cycle is the primary 
interest for engines, but this is a function of how well or badly all the energy exchanges over a 
cycle are performed. A more precise general understanding of these exchanges is the goal of this 
paper. 

As will be seen, the kinematic power ratios defined above are helpful in classifying and 
quantifying these energy exchanges for a given force process acting on a force-linear mechanism. 
Using these ratios, the theorem deduced below expresses the net work done by the y-actuator of a 
machine in terms of the work process applied to the x-actuator in such a way that the losses and 
gains can be clearly identified and easily understood. 

Let M = [. x,y,F(s,f,g)\ be a force-linear machine described by (3) and (f,g) a force process 
acting on it as 5 ranges over the closed bounded interval /, that is, / and g are functions of the 
state parameter s such that F(s,f(s),g(s)) = 0 for all s in I. It is assumed that all the functions 
here are piecewise continuous in s. Then the work output done by the y-actuator is 

K(g,y) = - / gdy = - J g(s)y'(s ) d s 

Let K be the union of all the closed subintervals of / in the interior of which x'(sj = 0, and let J be 
the closure of the compliment. Then J is a union of closed intervals and together with the intervals 
of K partition / with only endpoint overlapping. Thus 
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W 0 (g,y) = ~ / (g ~ go) dy - / g a dy = - / (g-g 0 )dy - / (g-g 0 )dy- / g 0 dy 
Ji Ji Jj Jk Ji 


( 5 ) 


The last term of (5) is a characteristic of the machine and does not depend upon the applied force 
process (f,g); this quantity will be denoted by W d \ 


w d = J gody 


The receding term in (5) arises from the portions of the working process when the x-actuator is 
stationary. This term represents a loss because when x'(s) = 0, relation (2) requires that 
(,g — g 0 )/(s) ^ 0. This loss depends upon the applied force process, but in practical machines it is 
usually a minor effect. In heat engines it represents work that must be supplied from the flywheel 
to carry the engine over periods where the piston is motionless. In most engines, this occurs only 
in two instants during each cycle, in which case this loss is zero. This loss will be denoted by W p : 


W p = I (g ~ go ) dy 
Jk 


Now the first integrand in (5) can be expanded by using the following calculation: 

~{g-go)y' = -[M+(f -fo) + -M-(f -foY)y' by(3) 

= [ k+ (f-f 0 ) + -k-if-fo)-}x l by(4) 

which holds when x'{s) fi 0. Now using the identity z = z + — z on x', k + , and l< , the last ex¬ 
pression expands to 


kl(f-f 0 ) + x' + - k~(f -fo) + x ,+ - k + _(f-fo)-x' + + k: (/ -/ 0 )-x' + + ktif-foYx' 

- kZ{f-f o yx'- - k + + (f-f 0 Yx'-+k-(f-f 0 Yx!- 


( 6 ) 


In this sum, the fourth and eighth terms are in fact zero under our assumption that s > 0. The 
fourth term 

kz(f-f 0 Yx' + 

could be positive only where 

k_ < 0, f — f 0 < 0, and x' > 0. 

if y 0, then the first and last inequalities imply A/ A :> 0, but this contradicts condition (2) when 
/ — f 0 < 0. If / < 0, then M_ < 0 and again (2) would not hold. Therefore the term is zero. A 
similar analysis shows the last term 

k-(f-f 0 ) + x'- 

is zero also for s > 0. 

In addition to assuming the functions x', /, M + and are piecewise continuous on the state 
parameter interval / = [a,b\, it is also reasonable to assume that they are of finite oscillation, that 
is, the inverse image of any open interval is a union of finitely many intervals open in I. This is a 
mathematical technicality and does not diminish the set of practical mechanisms to which the 
result can be applied. The assumption is sufficient for a mathematical proof that all the terms in 
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(6) are integrable; this proof is given in the Appendix. Therewith we have derived the following 
result. 

Theorem. If (/, g) is a force process acting on the x and y actuators respectively of a force-linear 
machine described by (3) with s > 0, then the net work done by the y-actuator is 

K(g,y ) = J Kif - f 0 ) + ±x + - ^ k~(f -f 0 ) + dx + - J ktif-foT d* + 

+ J k~(f ~ fo) dv - f kZ(f-f 0 )~ dx~ - J k+if-fot dx~ -W p -W d (7) 

Formula (7) yields insight into exactly where and when the desired work and the undesired losses 
occur in the operation of a machine. The first three terms refer to motion in the direction of 
increasing x since they contain dx + =x , (s) + ds. The first of these terms 

jK(f-fo) +dx+ 

represents a positive contribution to W 0 (g,y), the work output done through actuator y. These 
contributions occur exactly where 0 < k + ^ 1 and where the net input force f — f 0 agrees in sign 
with x'. Such processes are called efficacious. 

The next of the first three terms is 

- f k~(f - f a ) + dx + 

Because the integrand is non-negative and a minus sign precedes the term, it represents a loss 
relative to the desired work W 0 (g,y). These losses are incurred wherever k. is negative, / ^ f 0 , and 
x' > 0. Here the net force on the x-actuator is in the same direction as its movement, but the value 
of k + is negative and so work is required from the y-side to assist or carry the machine through 
such situations. This kind of loss occurs for example near the “dead centers” in a crank-con¬ 
necting rod mechanism. We will refer to this as a dead center loss although in general it can occur 
in other situations. 

The last loss term of this group is 

— J k-(f — fo) dx + 

A loss of this type occurs when the net x-side force, / —f 0 , opposes the motion of the x-actuator. 
We refer to such situations as forced processes since work is needed from the y-side in order to 
force the machine to move against the net input on the other side. Forced work processes occur in 
reciprocating engines whenever compression takes place above external pressure, or expansion is 
forced to take place below the external or buffer pressure [4], Note that when forced work occurs, 
the corresponding kinematic work ratio k is greater than unity. 

Interpretation of the three terms involving dr is similar. The term 

Jk+if-f 0 )- dx 

represents a gain in W 0 {g,y) through an efficacious process in which 0 < k d 1 ■ The term 
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Dead Center Process Efficacious Process Forced Process 


Fig. 9. Example of a force process f(s) acting on the x-side of a force-linear machine with x-preload f„. 


~ J k (/ - f° ) dx 

represents losses occurring when the net input force agrees with actuator velocity but k . < 0. This 
is another dead center type loss. Finally, the term 

-Jk+(f-f 0 ) + dx 

represents forced work; note k + ^ 1 when x' < 0. 

These losses and gains are easily identified graphically. Fig. 9 shows a force process / as a 
function of state parameter s for a force linear machine with a constant preload f„ on the x- 
actuator. The regions which correspond to positive contributions to W a (g,y ) are identified by the 
label G. These are the efficacious portions of the force process. Note that the work contribution to 
W a (g,y ) is not the area of the corresponding region, but is given by the first and fourth terms of 
(7). The regions corresponding to dead center type losses are identified by L D . The forced process 
loss regions are identified by L F . 

Because of the independence of k a from /, there is a general monotonic relation between the 
area of a region of one specific type between / and f a and the magnitude of the corresponding gain 
or loss term. To illustrate, just considering one region, say one of type L D , if one were to change 
the force process / so that the new region was contained (in the set-theoretic sense) within the 
original, then the L D loss for the new process would be smaller than for the original. This same 
type of monotonic relation holds for the other type regions as (7) makes clear. 


6. Application to reciprocating engines 

If an engine cycle with pressure p acts on a piston of area A attached to the x-actuator of the 
engine mechanism, and if p a denotes the external or buffer pressure on the non-workspace side of 
the piston, the force process applied is / = (p — p a )A and dx = x'(sjds = (\/A)dV. Thus the in¬ 
tegral terms in (7) have the “pdV” form 
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J I k «{ p ~ Po dVy wh ere ocJ,y <E {+,-} 

which directly relate to the pressure-volume diagram of the engine cycle. The first three terms of 
(7) in which y = + correspond to expansion processes. The three terms in which y = — relate to 
the compression processes of the cycle. 

Fig. 10 shows a simple pressure-volume diagram of a heat engine with a constant external or 
buffer pressure p a . Fig. 11 shows the corresponding force applied to the piston as a function of the 



Fig. 10. Example of a simple engine p-V diagram with constant buffer pressure. 



Fig. 11. Graph of piston force arising from the cycle of Fig. 10 vs. mechanism state parameter s. 



Fig. 12. The piston force process of Fig. 11 with the characteristics of a crank-connecting rod force-linear mechanism. 
In this example, there is no piston preloading. 
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Fig. 13. The p-V diagram of Fig. 10 with the mechanism characteristics of Fig. 12 superimposed. 

mechanism state parameter .v. Fig. 12 imposes on this the characteristics of a force-linear mech¬ 
anism typical of a crank-connecting rod system with no preloading. This diagram serves to 
identify how the engine cycle interacts with the engine mechanism. The portion of the cycle 
process from 1-2 produces a dead-center type loss. 2-3 is an efficacious process contributing work 
to the flywheel and output shaft. 3M is a forced expansion. 4-5 is more dead-center loss. 5-6 is an 
efficacious compression. 6-7 is another process incurring dead-center loss. Finally, 7-1 is a forced 
compression process. In simple cases such as this, most of the above information can be quali¬ 
tatively displayed on a slightly modified p-V diagram as shown in Fig. 13. 


7. Pressurization effects 

As an illustration of the utility of the above results, we consider here the effects of pressurizing 
or supercharging an engine. Most engine cycles, at least the textbook idealizations of them, are 
monomorphic [3], In a monomorphic engine, the workspace pressure function is proportional to 
the mass of the working fluid, everything else being the same. Examples of monomorphic cycles 
include the ideal Stirling, Carnot, and Otto with fixed temperature and volume extremes, more 
generally all similarly restricted Reitlinger and Crossley cycles, and many nameless others. 

The indicated work of a monomorphic cycle is proportional to the workspace fluid mass. This 
is the ideal response to charging an engine with more working fluid. However what counts in 
practice is not just the indicated output, but what comes out of the shaft. This does not always 
improve with charging. We treat here two important modes of charging: system charging, and 
workspace charging with buffering from below. 

7.1. System charging 

In engines having an enclosed pressure-worthy crankcase or buffer chamber, it is possible to 
pressurize both the workspace and the bufferspace uniformly as a system, by increasing the mass 
of the gas in each space by the same factor. Here it will be assumed that the bufferspace is 
monomorphic along with the workspace. A monomorphic workspace pressure p can be repre¬ 
sented mathematically by 
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p(s) = mpi(s) 

where m is the factor by which the mass of working fluid is increased in charging, and p\ (s) is the 
pressure function at a base charge corresponding to a unit factor, i.e., m = 1. Likewise the buffer 
pressure is representable by 

mpo 

where p a is the external pressure at base charge m = 1. Note that p„ need not be constant, but is 
usually nearly so in practice [5], 

The force applied to the piston at charge level m is therefore 

/ = m{pi - p 0 )A . 

Preloading on the piston, represented in our analysis by f a , is typically small in a well-designed 
and constructed engine. Taking f 0 = 0 as the best possible case, from (3) and (5) it follows that 

Wo(g,y) = - J^{g - go) dy -w d = m J (M_(p>\ -po)~ - M + (p x - p 0 ) + ) d y - W d . 

Recall that W d is a mechanism characteristic independent of the engine cycle and hence inde¬ 
pendent of m in the present application. This proves the following result. 

Theorem. The cyclic shaft work output of a monomorphic engine with a force-linear mechanism and 
negligible piston preload varies linearly with its system charge level. 

In practice, W d is relatively small, due only to preloads acting directly on the output shaft, such as 
contact friction from seals or constant speed viscous friction from lubricating or sealing fluid. In 
the limiting case when W d = 0, shaft output is directly proportional to system charge level m and 
mechanical efficiency is constant. 

7.2. Workspace charging with buffering from below 

The case where buffer pressure is left unchanged and the workspace charge level is elevated is a 
case of considerable practical importance since many engines are not equipped with a pressure¬ 
worthy crankcase or bufferspace. For such engines the buffer is the atmosphere and charging the 
workspace is the only possibility for improving performance through pressurization. 

Unlike system charging, workspace charging will not always result in an increase in cyclic shaft 
output. Once the charge level has been raised to the point where the entire cycle is above the buffer 
pressure level, shaft output will either ideally increase or decrease monotonically as the charge 
level is further increased. 

To see this, suppose m h is the charge level where the cycle is just buffered from below, so that 
mp\ > m h p\ ^ p 0 for all m > m b . Again regarding piston preload f as negligible, only the first, 
second, and sixth integral terms of (7) are possibly non-zero when m > m b because [mp x — p a ] = 0. 
Also taking W p = 0 = W d , shaft work is then 

W 0 {g,y) = J kf(mpi -p 0 ) + dV + - J kf(mp x - p 0 ) + dV + - J^ k+(mp x - p a ) + dV~ 
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These terms correspond to efficacious work on expansion, dead center loss on expansion, and 
forced work on compression respectively. Since [mp\ — p„] + = mp\ — p a when m > m h , shaft work 
can be written more simply as 


W 0 (g,y) 


k+(mp\ -p 0 )dV 


This shows that shaft output will vary linearly with charge level in the range [m b , oo). Whether 
Wo{g,y) increases or decreases as m is increased beyond m h is determined by the sign of the integral 


J k+pidV (8) 

Calculation of this integral provides a simple test of the response of an engine to charging above 
buffer pressure. Physically, the response depends upon the match between the mechanism char¬ 
acteristic k + and the base volume-pressure cycle (V.pi). A good match results in an engine in 
which workspace charging increases output; a poor match produces an engine in which super¬ 
charging degrades performance. Examples of engines exhibiting these behaviors are given below. 
In cases where the mechanism has simple constant effectiveness, much more can be deduced about 
an engine’s response to workspace charging [3]. 

As a simple example of the effects of workspace charging, consider an ideal Stirling cycle 
coupled to a force-linear mechanism having characteristics similar to a crank-connecting rod 
system. The engine is mathematically described by the following equations: 

*=\v(s) 

V(s) = V m + Vm 2 — (1 - coss) 


y = s 


M+(s) = (a - sin s) ( V M - V m )/{2A) 
M_(s) = {-a - sins)(E w - V m )/{2A) 
jiiRTh 

Pe = . for O^s^n 


Pc = 


V(s) 

mRTr 

nsi 


for n ^ s ^ 27i 


where 

s = crank angle, 

A = piston area, 

V (.v) = instantaneous engine workspace volume, 
V m = minimum workspace volume, 

Vm = maximum workspace volume, 

T h = maximum working fluid temperature, 

T c = minimum working fluid temperature, 
m = mass of working fluid, 

R = ideal gas constant, 



296 


J.R. Senft I Mechanism and Machine Theory 39 (2004 ) 281-298 


p E = pressure during isothermal expansion, p c = pressure during isothermal compression, 
a = mechanism constant. 

The mechanism in this example is mathematically simpler than the crank-connecting rod 
mechanism discussed above. Piston motion is sinusoidal. The constant a determines the effec¬ 
tiveness of the mechanism. The larger the value for a, the poorer the performance of the mech¬ 
anism; a = 0 corresponds to the frictionless case. Fig. 14(a) and (b) show the graphs of the 
kinematic power ratios for a = 0.1 and 0.3 respectively. It can be seen that the force characteristics 
of this model are quite similar to the actual crank-connecting rod mechanism. 

Fig. 15 shows the ideal Stirling cycle mathematically described above for the specific temper¬ 
ature and volume compression ratios 

T'c n <- , 

t = — = 0.5 and r = — =1.5 
Th V m 

and three working mass charge levels. At the lowest shown charge level, m =1.5, the minimum 
cycle pressure just equals the external buffer pressure depicted by the horizontal line. As the 
charge level is increased, the cyclic indicated work increases proportionally. With the buffer 
pressure held constant, shaft work will increase or decrease according to whether the sign of 
integral (8) is positive or negative. 

The base cycle corresponding to m = 1 has pressure function given by 



Fig. 14. Graphs of k + and k_ for a hypothetical force-linear mechanism with constant a = 0.1, and with constant 
a = 0.3. 



VOLUME VOLUME VOLUME 


Fig. 15. An ideal Stirling cycle with fixed temperature and volume extremes, constant external buffer pressure, and 
increasing working fluid mass. 
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The mechanism power ratio k + is given by (4) which in this example is 

k+ = (1 — U ] for s f 0,7i or and k+ = 0 at these points. 

\ sins/ 

Also, V'(s ) = Vu f Vm sins. Thus 


f*2n 


k+p\dV = / k + (s)p\(s)V'(s) ds = 


(r — l)RT} 


H 


2(1—t)lnr U7r(l+T) 


1 


V~r 


The factor in the brackets above works out to 0.426 for the chosen cycle with t = 0.5 and 
r = 1.5 with the better mechanism having a = 0.1. This positive value means that shaft output will 
increase as workspace charge is increased. For the poorer mechanism having a = 0.3, the factor 
calculates to -0.343 which means that engine output will steadily decrease as the workspace 
charge is increased above buffer pressure. 


8. Conclusion 

The force-linear concept introduces a simple yet realistic model of machines appropriate for 
analysis under conditions of fixed or constrained motion. Applied to reciprocating engines, it 
permits the classification and quantification of the various frictional losses suffered by an engine in 
a way which improves intuition on how an engine cycle and a mechanism interact to determine 
overall engine performance. The force-linear model is also useful for determining the performance 
response of an engine to changes in its thermodynamic cycle, such as those due to supercharging. 
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Appendix 

Lemma. The functions lfx' y are piecewise continuous for all choices of a,/f,y £ {+, —}. 

Proof. From the assumptions of piecewise continuity and finite oscillation, it is possible to find a 
partition P : a = s 0 < si < s 2 < ■ ■ ■ < s„ = b of 1 into interior-disjoint subintervals f = [.$;,■ i. .v,] 
with the following properties: 

• the functions x\ /, M + and M are continuous on the interior of each f and are continuously 
extendable to f, 

• one and only one of the following holds throughout the interior of each If. 

x' >0, x 1 = 0, or x' < 0 

• k + and k are continuous on the interior of each /,, 
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• for each a e {+, —}, throughout the interior of each Ij, exactly one of the following holds: 
k x < 0, O^A^l, or 1 < k a . 

Using this partition, it is easily seen that all the terms in expression (6) are piecewise continuous 
functions of .v. On any of the subintervals /, of the partition P, consider the function k 7 (s)x'(sy 
where a, y e {+, —}. In a case where x'(s) is identically zero on Int(/ y ) so also is k x (s)x'(s)‘ zero on 
Int(/ ; ). In the other cases, jc'(s) is of one sign throughout the interior of Ij . Therefore x'(s) y is jc'(s), 
0, or — x'(s) throughout Int(/ ; ). Hence x ly /x' is identically equal to one of the constant functions 1, 
0, or -1 in Int(/ ; j; denote this constant by c. Then 

k a (s)xf(s) y = -M a (s)y' (s)^jrl- = -M x (s)/(s)c 

x'(s) 

in lnt(7 ; ). The right hand side of this equation is continuous in lnt(7 ; ) and the limits from the 
interior to the endpoints exist because they do for each factor. This proves that the functions 
k a (s)x'(s) y are piecewise continuous on I for all choices of a,y 6 {+,—}• From this the lemma 
follows because 

[k a x' y } + = k^[x' y } + + k-[x' y Y = k^x' y 

and 

mt= k w+ k r ] + = D 
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